Nonequilibrium adsorption of 2AnB patchy colloids on substrates 
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We study the irreversible adsorption of spherical 2AnB patchy colloids (with two A-patches on the 
poles and n _B-patches along the equator) on a substrate. In particular, we consider dissimilar AA, 
AB, and BB binding probabilities. We characterize the patch-colloid network and its dependence on 
n and on the binding probabilities. Two growth regimes are identified with different density profiles 
and we calculate a growth mode diagram as a function of the colloid parameters. We also find 
that, close to the substrate, the density of the network, which depends on the colloid parameters, is 
characterized by a depletion zone. 



I. INTRODUCTION 

Patchy colloids have functionaUzed surfaces yield- 
ing directional colloid-colloid interactions. [THS] This fea- 
ture opens the possibility of controlled bottom-up self- 
organization of novel materials [SUlOj and it has moti- 
vated the development of several experimental techniques 
to produce patchy colloids. [TTHT5] Theoretical studies 
have been focused on the thermodynamic properties with 
several numerical techniques being adapted to tackle the 
patch-colloid cooperative behavior in solution. [21 IT9H2T] 
A range of models for the colloids have been extensively 
studied and used to represent many complex systems, 
such as, amphiphilic molecules, colloidal clays, proteins, 
and DNA nano-assemblies.[2l El E H EQI [^MS] These 
studies have revealed that the control over the valence 
leads to rich phase diagrams with new self-organized 
phases and interesting properties, such as, glass transi- 
tions, colloidal gels, and bigels. p7H5^ 

Studies of equilibrium properties provide insight on the 
possible thermodynamic structures and their stability. 
However, experimentally, equilibrium structures are very 
often not accessible. In particular, in the limit of irre- 
versible binding (considered here) , the network of colloids 
is determined by the time sequence of patch-patch bond 
formation. |33l |34| In general, these kinetically trapped 
structures are different from the thermodynamic ones. 
The interest on the identification and characterization 
of these structures is twofold: designing experimental 
strategies to avoid them when the equilibrium structures 
are the goal and the search for novel material properties. 

The presence of a substrate improves the degree of con- 
trollability over self-organization of patchy particles even 
under equilibrium conditions. '35( l36] Recently, the study 
of the nonequilibrium adsorption on a substrate of two- 
and three-patch colloids revealed the existence of an op- 
timal fraction of two-patch colloids where the density of 
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the film is maximized. [33j This result is in contrast to 
the equilibrium structures where, a monotonic decrease 
of the density is observed. [22] Several questions about 
nonequilibrium patchy-colloidal systems are still open, 
namely, the effect of confinement [3], the inclusion of 
patterns on the substrate, [37l [38] or the possibility of dis- 
tinct interaction between patches. |39) In this work, this 
last question will be investigated in detail. 

A particularly interesting model considers patchy col- 
loids with distinct patch-patch interactions, namely 
2AnB (two patches of type A and n patches of type 
B), also known as Lisbon colloids. Even the simplest 
version of this model {n — 1 with all patches evenly 
distributed) yields interesting properties in the bulk, 
such as vapour-liquid transitions and a vanishing crit- 
ical point. The liquid- vapour critical point was stud- 
ied as a function of the AA, AB, and BB interaction 
energies. [401 141] and a rich number of network structures 
was found to be responsible for this unexpected phase 
diagram. Both numerical and theoretical results confirm 
the unusual thermodynamic and percolation properties 
of this model. ^J2l |33] More complex 2AnB colloids with 
two patches of type A located on the poles and nine 
patches of type B located equidistantly on the equator 
have been studied. [39|, 144] An equilibrium liquid- vapour 
binodal was found to be re-entrant with the liquid (a 
network fluid with a large number of branches) and the 
vapour (a fluid of linear chains with little or no branch- 
ing) densities vanishing as the temperature approaches 
zero. This is in line with the phase diagram predicted by 
Tlusty and Safran for dipolar fluids. [31] This re-entrant 
behavior was also found on a lattice model in two and 
three dimensions. [461 H71 

Until now, the study of 2AnB colloids has been re- 
stricted to equilibrium properties. This gives an idea of 
the numerous possibilities that can be expected in the 
richer nonequilibrium systems. Following this lead, in 
this work, we will study the nonequilibrium adsorption 
on an attractive substrate, using our model, proposed 
in Ref. [33j. where the bonding between colloids is irre- 
versible. Under this constraint, the system never reaches 
an equilibrium state. As in previous models of 2AnB 



colloids, [Sni m] we consider two patches A in the poles 
and n patches B around the equator. An example for 
possible configurations is shown in Fig. [ija) with two 
schematic representations of 2AnB colloids with n — 4 
and n = 9, respectively. We study the effect of varying 
both the number of patches of type B and the interaction 
between patches. We find a promising control over the 
film structure by tunning the 2AnB colloids properties. 
A detailed analysis of the film properties is presented, 
based on kinetic Monte Carlo simulations. We observe a 
wide range of coverages of the first layer by changing both 
the number and the interaction between patches, and the 
emergence of a density depletion zone near the substrate. 
A new growth mode diagram depending on the number 
of i?-patches and sticking coefficients is proposed. 

This article is organized as follows: In Sec. [ll] we in- 
troduce the model and the underlying physical considera- 
tions, as well as the simulation parameters; In Sec. |III| the 
results for the region close to the substrate are presented 
and the density depletion in this region is discussed; In 
Sec. |IV| the film growth modes are discussed; Finally, in 
Sec. I VI we draw some conclusions. 



II. MODEL 

We consider the stochastic model for the adsorption 
of patchy particles on a substrate, recently introduced 
in Ref. [33] . Patchy colloids are typically described as 
spherical particles with patches distributed on their sur- 
face, with both a short-range repulsive interaction and 
a patch-patch attractive interaction. We assume chemi- 
cal bonding between patches (highly directional and irre- 
versible within the timescale of interest) and we describe 
the patch-patch binding as a stochastic process, where a 
bond is established in two different situations: when the 
patch of the adsorbing colloid is aligned with the patch 
of a previously adsorbed colloid or when the patch-patch 
interaction promotes the alignment of the landing col- 
loid with the available binding site. An interaction range 
around the patches is defined such that two colloids ef- 
fectively interact (and eventually establish a bond) when 
their interaction ranges overlap. Here, we fix the interac- 
tion range within an angle lower than 7r/6 rad relative to 
the center of the patch. [JS] This sets a maximum interac- 
tion TTa/6 between patches on the surface of the colloids, 
where a is the diameter of the colloid, which we consider 
unity without loss of generality. 

To describe the colloids motion in solution, a Brown- 
ian dynamics is considered where the collision with the 
solvent is taken as a Poisson process, such that the time 
between collisions. At, follows an exponential distribu- 
tion. 



p{At) = exp(-i?At), 



(1) 



where R is the collision rate. At each colloid/solvent col- 
lision a new velocity is assigned with the direction uni- 
formly distributed and with absolute value v following 
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FIG. 1. Schematic representation of 2AnB colloids, where A 
patches are in red and B patches in green, (a) With n — 4 
(left) and n = 9 (right), (b) Only AA (left) or AB (right) 
bonds can be established, i.e., BB bonds are neglected, (c) 
Representative structures formed by assembly, with chains of 
both AA and AB bonds (left most structures), and junctions 
of bonds, i.e., multiple bonding of the same colloid (more 
than two bonds). Junctions are of type {AA)x(AB)y, where 
a; = 0, Na and y = 1, Nb + Na {x + y > 3). 



the Maxwell-Boltzmann distribution at the thermostat 
temperature, given by. 



p{v) 



1 



2ttT 



exp 



2T j 



(2) 



with T the thermostat temperature in units of kg/m, 
where m is the mass of the colloid and kg the Boltzmann 
constant. With this model, it is possible to fine tune the 
diffusion coefficient D, with the proper combination of T 
and R. For the sake of simplicity, we consider a constant 
diffusion coefficient of O.lcr^/s. 

In the presence of a substrate, there are two character- 
istic time scales: the inter-arrival time of colloids at the 
substrate (typically related with the fiux); and the bind- 
ing time, defined as the time necessary for the chemical 
bond to be established. At low concentration of colloids, 
the inter-arrival time is considered much larger than the 
binding time. Therefore, we can assume that colloids ar- 
rive one at a time to the substrate and adhere instanta- 
neously. For the purpose of this study, the interaction of 
a colloid with the substrate is isotropic and they adsorb 
with no defined direction. 

The process of establishing a chemical bond is a ther- 
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mally activated process, characterized by an activation 
barrier. In general, this barrier is different for AA, AB, 
and BB bonds. Here, following previous studies of 2AnB 
colloids, [Snilll] we consider that the activation barrier for 
BB is so much higher than the others, that such bonds 
can, in practice, be neglected. We define the binding 
probability for AA and AB bonds as P/i^ and Pab- For 
simplicity, we consider the limit Paa = 1 and define the 
sticking coefficient tab — Pab /Paa- The two possible 
bonds are represented in Fig.[l]jb). 

The parameters n and tab are related. The larger 
their values, the higher the probability that an AB bond 
is formed. Yet, as we show below, a different quanti- 
tative and qualitative dependence on each parameter is 
observed. Figure [ijc) shows some possible local arrange- 
ments obtained with 2AnB colloids. The examples are: 
a sequence of consecutive AA bonds, always resulting 
in a linear chain of colloids; yli?-chains, formed by se- 
quences of AB-bonds, which are rarely linear (linear AB- 
chains are only geometrically possible for even values of 
n) whose shape depends strongly on n; and {AA)x{AB)y- 
junctions, where x — 0,...,Na and y 1,...,Nb + Na 
(x + y > 3), which are colloids with more than two bonds. 

In a previous work for colloids with identical 
patches, [33] we found three distinct regimes in the den- 
sity profile: the surface layer, where the effect of the 
substrate, on the colloid network, plays a major role; 
the liquid film, where the density of the film reaches a 
constant value during growth; and the intcrfacial region, 
which is mainly related to the roughness of the film. In 
this article, we focus on the surface layer and liquid film 
for 2AnB colloids. 



III. SURFACE LAYER 

The layer of colloids directly in contact with the sub- 
strate is the precursor film. As we show here, for a 
patchy-colloid network with strongly anisotropic colloid- 
colloid interaction, the organization of the first layer 
affects the growth of the remaining structure. For 
nonequilibrium adsorption of patchy colloids with identi- 
cal patches, it has been shown [55] that the density of the 
first layer remains constant, independently of the system 
size, but increases with the diffusion coefficient D of the 
colloid. 

The coverage of the first layer, i.e., the number of par- 
ticles per unit of area, as a function of n and tab is 
shown in Figs. [2]Ja) and (b), respectively. It is observed, 
in Fig. [2ja), that for tab ~ 1, the coverage of the first 
layer decreases with n, while for tab ^ 1, it remains 
constant. At low tab the probability of binding to a pre- 
viously adsorbed particle is very low and, consequently, 
particles are able to diffuse towards the substrate and 
bind directly to it. This binding is solely limited by the 
colloid-colloid excluded volume interaction and, in the 
limit of vanishing tab i the coverage of the substrate by 
particles directly in contact is expected to resemble the 
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FIG. 2. (a) Coverage of the substrate, 9 (number of particles 
directly adsorbed on the substrate per unit area), as a func- 
tion of the number of _B-patches, n, for values of tab between 
O.Of and 1. (b) Coverage of the substrate, 6, as a function of 
the sticking coefficient, tab, for values of n = {f , 2, 3, 4, 5, 6}. 
(bottom) Snapshots of the first adsorbed layer (directly in 
contact with the substrate) for three different points on the 
plots (a) and (b), namely, f. A, and F/6, where the depen- 
dence of the density on n and tab is observed. Simulations 
performed on a substrate of lateral size 32 in units of particle 
diameter. 



Random Sequential Adsorption model. [49] with a jam- 
ming coverage, ^oo = 0.6969. The effect of the number 
of patches will then be negligible for tab ^ f (like, e.g., 
TAB = 0-01 in Fig. [fa)). 

In general, for any value of tab, a decrease of the cov- 
erage of the first layer with n is observed. This is related 
with the increase in the probability that a particle binds 
to the network, thus enhancing the screening effect. A 
similar qualitative behavior is observed with increasing 
rAB ■ This dependence on n and tab is also visible in the 
snapshots of the first layer, Fig. [2| for limiting values of 
n and tab , where the density at low n and tab is higher 
than at higher values. 

For patchy colloids with identical patches, a decrease 
of the density with the distance to the substrate is ob- 
served for the surface layer. [33] However, this is not the 
case for 2AnB colloids, with distinct patches. Figurejsj^a) 
shows several density profiles obtained with these parti- 
cles, where a minimum is observed for z < 2. The pres- 
ence and position of this minimum is robust with respect 
to the values of n and tab- Below, we discuss this be- 
havior in detail. 

Figure |3][a) shows the density profile for some values of 
rAB and n where the minimum is observed. An explana- 
tion may be drawn from Fig. [sj^b) , where a density profile 
is shown, with a binning of Az = O.f . Here, a narrow 
peak appears at a distance of f .3 < z < f .4 from the sub- 
strate (z/ in the figure). This is a direct consequence of 
the high density of the first layer and the low probability 
of adsorption on _B-patches. The scheme A in the bottom 
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FIG. 3. (a) Density profile of the colloidal network, close 
to the substrate, from 2: = to z = 10 and various values 
of n and tab, with a binning of Az = 1. From bottom to 
top {Nb = I, TAB = 0.01}, {iVfl = l,rAs = 0.1}, {TVs = 
6, Tab = 0.1}, and {Nb = 6,rAB = 0.01}. (b) Density profile 
of the colloidal network on a substrate, with binning of Az = 
0.1 from z = to 2 = 4, for the same parameters and symbols 
of plot (a). The shadowed region corresponds to the volume 
used for the calculation of the density at 1 < 2 < 2 in plot 
(a). Labels A, B, and C indicate the first three peaks in the 
density above the layer directly in contact with the substrate, 
(c) Density of the region at 1 < 2 < 2 as a function of vab 
for n — 1 (open squares), 2 (solid squares), 3 (open circles), 
4 (solid circles), 5 (open triangles), and 6 (solid triangles). 
Simulations performed on a substrate of lateral size L = 32, 
in units of particle diameter, and a total of 61440 deposited 
colloids, (bottom) Schematic representation of the structures 
at the peaks marked in (b). 



of Fig.|3]illustrates the most common configuration of the 
colloid connected to the one that is directly bound to the 
substrate. For an isotropic orientation of colloids in the 
first layer, the probability of finding a colloid near the 
normal to the substrate is very low and, due to excluded 
volume, adsorption of colloids on patches that form an 
angle lower than 7r/2 rad with the substrate is also lim- 
ited. Considering a coverage around 0.65 particles per 
unit of area, we obtain an area of approximately 1.54cr^ 
per particle. The average distance between particles can 
be estimated from the diameter of such regions, yield- 
ing Ar « 1.4(7. The minimum height of a colloid, not 
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FIG. 4. (a) Density difference, p*(0) — p*(l), correspond- 
ing to < 2 < 1 and 1 < z < 2, respectively, as a func- 
tion of Tab for values of n = {1,2,3,4,5,6}. (b) Density 
difference, p*(2) — p*(l), corresponding to 2 < 2 < 3 and 
1 < 2 < 2, respectively, as a function of tab for values of 
n — {1, 2, 3, 4, 5, 6}. (bottom) Snapshots for, from left to right 
and top to bottom, {Nb = 3, tab = 0.01}, {A^'b = 6, tab = 
0.01}, {iVB = 3, TAB = 1}, and {Nb = 6, tab = 1}. Results 
averaged over 100 samples on a square substrate of lateral size 
L = 32 in units of the particle diameter and a total of 61440 
deposited colloids. 



directly in contact with the substrate, would correspond 
to a particle resting on top of two other particles at this 
distance. In that case, a minimum z 1.21 is obtained. 
Since bonds can only be established along the direction 
of the patches, the maximum is observed in the range 
1.3 < z < 1.4. 

What enhances this behavior with 2AnB colloids is the 
lower probability to bind to B-patches, since BB bonds 
are absent and tab < 1- However, when an attempt 
of adsorption on a i3-patch is successful, the new bond 
works as the seed to a lateral chain, which consequently 
screens the adsorption on the region 1 < z < 2 (scheme 
B in the bottom of Fig. |3| . This few successful adsorp- 
tions on lateral B patches are represented in Fig. [sf^b) 
as a small peak just bellow z = 2 (zu). When r^B w 1 
the colloidal network resembles the one for colloids with 
identical patches, where the density decreases monotoni- 
cally with the distance to the substrate due to branching. 
In this limit, the effect of the dense coverage of the first 
layer is negligible, as new colloids can adsorb equally on 
any patch. 



Figure 
in Fig.^ 



STc) shows the density of the shadowed region 
d) as a function of n, for different values of rj^g. 



For lower tab, the density increases with n. By contrast. 



for higher tab, it decreases with n. For low values of tab 
the increase in the density is related to larger branching 
and lateral growth of the film, since the first layer has 
a constant density regardless of n, as stated above. The 
higher the value of n the higher the branching rate and, 
consequently, the number of lateral adsorptions on col- 
loids of the second layer. The number of possible bonds 
along the equator is limited by the excluded volume in- 
teraction among colloids. Thus, it is expected that above 
n = 4, the efficiency of adsorption does not improve with 
n. For high values of tab, the density of the shadowed 
region in Fig. [sj^b) varies with the density of the first 
layer. 

The control of the density in the depletion zone can be 
achieved through the choice of n and tab (see Fig. |4]). 
The density close to the substrate (z < 1) is always larger 
than in the second layer (1 < z < 2). Due to the higher 
probability to stick to the substrate than to any patch 
(see Fig. |4|^a)). The second layer (1 < z < 2) can have 
a density higher or lower than the third (2 < z < 3) de- 
pending on n and rAB- This is a direct consequence of 
the distinct patch-patch interaction since, as observed in 
Fig. |4]jb), for higher values of tab, the depletion dis- 
appears, and the behavior is similar to that observed 
for colloids with identical patches. [33] Snapshots of the 
film close to the substrate show a reduced density for 
1 < z < 2, more evident for low values of tab (see 
top snapshots in Fig. |4]) , and the effect of lateral chains 
(scheme B in Fig.|3|. 



IV. LIQUID FILM 

For the adsorption of colloids with identical patches, it 
was observed that the lateral growth of the colloidal net- 
work leads to a saturation in the film density (p = pi), 
due to the finite size of the substrate (see right snap- 
shot in Fig. [5]). Surprisingly, with distinct patches, we 
found for low values of rAB a new type of growth where 
no saturation of the density is observed (left snapshot in 
Fig. 5]). Figure [5] shows snapshots and density profiles 
for the two different growth modes, where the difference 
in the film structure is visible. For higher values of tab 
we observe a higher branching rate due to the adsorption 
on both A and B patches (bottom zoom ofFig.[5|. The 
junctions consist mainly of v4_B-bonds, which can initi- 
ate the formation of ^i3-chains (Fig.jljc)). These chains 
might have different geometries as pictured in Fig. [ijc). 
By contrast, for lower values of tab, the ratio and size of 
AA-chains (Fig. Iljc)) increases with the number of de- 
posited particles (top zoom of Fig.js]), which is associated 
with a decrease in the density, since a compaction of the 
film is hindered by the linear AA-chains. 

The different behavior at tab ^ 1 is a consequence of 
an increase in the relative density of A A over AB bonds. 
By contrast with 2AnB colloids in equilibrium, where 
bonds are reversible and a low value of tab leads to a 
negligible number of j4i?-bonds, under nonequilibrium 
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FIG. 5. Density profile of the colloidal network on a substrate 
showing three different regimes: surface layer (green, bot- 
tom), liquid film (red, middle), and interfacial region (blue, 
top). Left: Reduced density, p*{z) (number of particles per 
unit of volume), as a function of z, for n — 1 and vab ~ 0.01; 
Right: Density, p*{z) as a function of z, for n = 1 and 
rAB ~ 1; (center) Magnified region of the colloidal network 
(top) n = 1 and tab ~ 0.01, and (bottom) n = 1 and 
rAB = 1. Color gradient indicates AA-chains (red) and AB- 
chains (green), with junctions represented as black spheres. 
Both cases after the adsorption of 40960 particles, averaged 
over 100 independent realizations. Zmax/2 is the value of half 
the maximum thickness of the film. 



conditions, even for low rABi A_B-bonds still contribute 
to the structure of the film. In equilibrium, for low values 
of rAB, the structure consists mostly of AA-chains[39] 
at an almost constant density around the equilibrium 
value. Out of equilibrium, for a network of 2AnB col- 
loids in the limit of rAB ^ 1, the density of bonds of 
type AB decreases with the number of adsorbed parti- 
cles. This effect is due to the increasing dominance of AA 
chains, which will screen available patches inside the film 
and, consequently, keep the density of the film low. The 
large number of B patches available for an A_B-bond are 
shadowed by particles above, and remain unreachable for 
newly incoming colloids. 

For colloids with identical patches, the liquid regime 
was defined as the regime, above the surface layer, where 
the density saturates. [33] Here, we found a growth mode 
where the density of the film does not saturate. Despite 
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a continuous decrease of the density (left of Fig. [s]) , a liq- 
uid regime can be identified as the region where the finite 
size of the substrate plays a relevant role. This decrease 
of the film density with the number of deposited particles 
is a result of the dominance of the rate of formation of 
v4A-chains over the rate of formation of junctions or AB- 
bonds. This behavior is evident from the derivative of 
the density of AA- and AB-honds, 



dz 



and 



dNA 



— . We name this growth 



in Fig.Je^a), where > 
mode Chain Dominated Growth (CDG). This is in sharp 
contrast with the second growth mode, where the density 
saturates (right of Fig. ^ and = = 0. We 

name this mode Junction Dominated Growth (JDG). To 
compute the derivative of the density of AA- and AB- 
honds, the liquid regime region needs to be defined pre- 
cisely. For the JDG mode, the liquid regime is the re- 
gion where the density of bonds remains constant. In 
the CDG mode, as a result of the decreasing density, 
the density of bonds exhibits a maximum (see Fig.|6ja)). 
As we go away from the substrate, the density of bonds 
initially increases due to the absence of binding to the 
substrate and the low density of the depletion zone close 
to the substrate. The maximum can then be considered a 
transition point between the surface layer and the liquid 
regime, and can be used as a reference for the measure- 
ments in the CDG mode. 

To characterize the growth-mode diagram and its de- 
pendence on n and tab we compute the derivative of 
the density of A A and AB bonds as a function of z. In 
Fig. one can see that for the CDG mode the den- 
sity of AB bonds decreases with z, and for lower values 
of n the density of AA bonds also decreases. The over- 
all decay of both can be related to the decrease of the 
film density. In the CDG, however, the rate of decrease 
of ^-B-bonds is higher than that of vlA-bonds. We can 
then characterize the growth modes as, 

dNAB ^ dNAJ 



dz 



and 



dNAi 



CDG 



JDG. 



(3) 



dz dz 

A growth diagram is shown in Fig. [6]jb) for the range 
of n and tab investigated. For tab < 0.1 the growth 
is CDG and for tab > 0.1 it is JDG. For n = 1 and 
fAB = 0.1 the growth mode is still CDG. 

For n < 4, rAB has a large effect on the formation 
of junctions, however this effect disappears for n larger 
than four due to steric limitations. In the limiting case 
of Tab li there is a strong dependence of the overall 
density on n, so, for lower values of n, the number of 
unsuccessful binding attempts is higher and the density 
increases through the diffusion of colloids into the fjords 
of the film. For the case of tab ~ 1, a higher number 
of successful bonds leads to a larger ramification of the 
film. In the limit of rAB ~ 1, n will have an effect similar 
to rAB, i-e., it increases the rate of bonding, as most 
of the patches are not available for binding due to the 
screening by newly deposited colloids. Since for CDG no 
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FIG. 6. Analysis of the growth modes for different values 
of n and rAB- (a) Number of AA- (red open squares) and 
yl-B-bonds (green solid squares) per unit volume. Neon * {z), 
as a function of z, for values of tab and n from left to 
right: {tab = 0.01; A^B = 1}, {tab = 0.01; A^b = 6}, 
{rAB = 1;A''b = 1}, and {rAB = 1',Nb = 6}. Solid lines in- 
dicate the derivative of the function in the liquid film regime, 
(b) Growth mode diagram for values of n as a function of 
rAB, for the Chain Dominated Regime (red open squares, left 
blue region) and the Junction Dominated Regime (green solid 
squares, right orange region). Results for 61440 deposited 
particles on a substrate with L = 32 averaged over 100 sam- 
ples. 



saturation of the density is observed, the density at half 
the thickness of the film, p(zma2;/2) was measured (see 
Fig. [?]). Taking z^axj'^ as defined in Fig. [5] for n > 2, 
the density decreases with tab ■ This is a consequence of 
a higher ratio of A_B-bonds, which favors the formation 
of junctions and A_B-chains, and, consequently, promotes 
ramification. Thus, it results in the appearance of regions 
of very low density. For n < 2, the density exhibits a 
maximum at the transition between the CDG and JDG 
modes. This is due to the effect of AA-chains, which, 
for the CDG mode, dominates at lower n. We note that. 
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FIG. 7. Density at z = Zmaxj'i as a function of tab for 
different values of n: one (open squares), two (solid squares), 
three (open circles), four (solid circles), five (open triangles), 
and six (solid triangles). Results for 61440 deposited particles 
on a substrate with L = 32 averaged over 100 samples. 



in the CDG mode, the density of the film is expected to 
vanish in the limit z — >■ cxi. 

The results of Fig. [Tj at values of r^s, inside the JDG 
mode, reveal a non-uniform dependence of the density 
on n. This unusual behavior is observed at values of 
Tab = 0.1 and tab = 0.2. Due to large fluctuations of 
the density at Zmax/"^ (from the error bars of Fig. [t]), 
the saturation density pi was chosen as the variable to 
characterize the film. Figure [sj^a) shows that the film 
density exhibits a maximum that depends on tab and n. 
For TAB =0.1 the maximum is found at n = 4 and for 
Tab = 0.2 at n = 2. For these pairs of tab and n, the 
film grows according to the CDG mode. In this growth 
mode, favored by junctions, the maximum in the density 
is driven by the maximum number of bonds per parti- 
cle. This is the case for n = 4, as seen in Fig. 8]^b), where 
the average number of bonds per particle, with more than 
one bond, is plotted. We observe that for both tab =0.1 
and TAB — 0.2 the maximum in the average number of 
bonds occurs for n « 4, implying that the configura- 
tion with two patches in the poles and four around the 
equator (octahedral geometry) maximizes the formation 
of junctions. As a conclusion, an increase in n does not 
necessarily promote junctions, which explains the max- 
imum in Fig. [sjja). For tab — 0.2, however, branching 
and formation of voids shift the density maximum to- 
wards lower values of n (around two). 



V. CONCLUSIONS 

We have studied the structural properties of colloidal 
networks obtained from the irreversible adsorption of 
2AnB colloids on a substrate. This is the first study 
of nonequilibrium growth of patchy colloids with a het- 
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FIG. 8. Analysis of the liquid film density near the transition 
between growth regimes, (a) Liquid film density, pi, as a. func- 
tion of n for Tab = 0.1 (open squares) and tab = 0.2 (open 
circles), (b) Average number of bonds per colloid (colloids 
with more than one bond), < Ni,ond >, as a function of n for 
Tab =0.1 (open squares) and tab = 0.2 (open circles). Re- 
sults for 61440 deposited particles on a substrate with L = 32 
averaged over 100 samples. 



erogeneous distribution of patches, namely, two types of 
patches restricted to the poles and the equator. We have 
shown that these 2AnB colloids are endowed with new 
tunable parameters for the control of the network prop- 
erties. 

We found that 2AnB colloids exhibit different struc- 
tures from those observed in (2 -I- n)A colloids both at 
the surface and within the liquid film. A depletion of 
particles at one particle diameter from the substrate is 
observed. We also show that the coverage of the substrate 
can be controlled through the number of patches in the 
equator and the unlike patch-patch interaction. For sim- 
plicity, we have considered an isotropic particle-substrate 
interaction. In general, some anisotropy is expected due 
to the presence of the patches. We have shown[33] that, 
when only one type of patch is considered, the liquid 
film regime is not affected by the details of the interac- 
tion with the substrate. Nevertheless, differences have 
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been observed in the surface-layer regime. In the pres- 
ence of two types of patches a richer surface-layer regime 
is expected depending on the selectivity of the particle- 
substrate interactions. 

We summarized the different growth modes in a two- 
parameter diagram, namely n and tab- We have identi- 
fied two growth modes that we named Chain Dominated 
Growth (CDG) and Junction Dominated Growth (JDG). 
The structure of the film in the two modes is intrinsically 
different, for example: in the CDG mode the density of 
the film decreases with the thickness while for the JDG 
mode it remains constant. We also propose to control the 
overall film density by tunning the number of patches and 
their interaction. 

These networks might be of relevance in the fields of 
microfluidics and filtering. For this purpose, the exper- 
imental realization of 2AnB colloids is necessary and 
could be achieved by using triblock colloids 8 or colloidal 
particles functionalized by DNA. 17, 50 As in previous 
work,[33_ we have focused on the nonequilibrium prop- 



erties in the limit of irreversible binding, usually associ- 
ated with low temperatures, where the obtained colloidal 
structures do not change over time. However, on very 
long timescales or at higher temperatures, the possibility 
of breaking bonds can no longer be neglected and the film 
is expected to relax towards the equilibrium structure. 
The time evolution of these structures and the identifi- 
cation of the characteristic timescales of relaxation are 
interesting open questions which are left for future work. 
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